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I ' Abstract 
(D , 

(— I , A model for the Pomeron at t = is suggested. It is based on the idea of a 

finite sum of ladder diagrams in QCD. Accordingly, the number of s-channel 

^ , gluon rungs and correspondingly the powers of logarithms in the forward 



scattering amplitude depends on the phase space (energy) available, i.e. as 
energy increases, progressively new prongs with additional gluon rungs in 
the s-channel open. Explicit expressions for the total cross section involving 
two and three rungs or, alternatively, three and four prongs (with In^(s) and 
In'^(s) as highest terms) is fitted to the proton-proton and proton-antiproton 
total cross section data in the accelerator region. 
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1 Introduction 



It is widely accepted that the Pomeron in QCD corresponds to an infinite gluon 
ladder with Reggeized gluons on the vertical lines (see Fig. 1), resulting [|T], ^ Q 
in the so-called supercritical behavior at ~ s"'^^\ where a(0) is the intercept of the 
Pomeron trajectory. However, at finite energies only a finite number of diagrams 
contributes. The lowest order diagram is that of two-gluon exchange, first consid- 
ered by Low and Nussinov 0|. The next order, involving an s-channel gluon rung, 
was studied e.g. in the papers [@, ^. The problem of calculating these diagrams is 
twofold. One problem is connected with the non-perturbative contributions to the 
scattering amplitude in the "soft" region. It may be ignored by "freezing" the run- 
ning coupling constant at some fixed value of the momentum transfer and assuming 
that the forward amplitude can be cast by a smooth interpolation to t = 0. More 
consistently, one introduces a non-perturbative model of the gluon propagator 
valid also in the forward direction. The second problem is more technical: at any 
given perturbative order a", the leading contribution in the s ^ oo limit, propor- 
tional to (Q;sln(s))", is given by a subset of all the Feynman diagrams contributing 
at that perturbative order; each of these diagrams consists of a leading term in the 
s — s> oo limit and of a non- leading, negligible part. The leading contributions from 
all orders in perturbation theory can be resummed ^ ^. For non-asymptotic en- 
ergies, however, at any order in the coupling constant subleading terms are present 
coming both from the neglected diagrams and from the neglected part of the leading 
diagrams. Although functionally the result is always the sum of increasing powers 
of logarithms, the numerical values of the coefficients entering the sum is lost unless 
all diagrams are calculated. 

The summation and convergence of an infinite series is a known problem in 
physics. As discussed in a recent paper 0, various situations may occur, where a 



finite series approximates the exact result better than the infinite sum does. Since, 
as stressed above, the coefficients of the perturbative series are not known from 
QCD even for t = (their calculation in the non-forward direction is much more 
tricky), the convergence of the series is also unknown. 

Conversely, one can expand the "supercritical" Pomeron ~ s"^^^ in powers of 
ln(s). Such an expansion is legitimate within the range of active accelerators, i.e. 
near and below the TeV energy region, where fits to total cross sections by a power 
or logarithms are known |^ to be equivalent numerically. Moreover, forward scat- 
tering data (total cross sections and the ratio of the real to the imaginary part of 
the forward scattering amplitude) do not discriminate even between a single and 
quadratic fit in ln(s) to the data. 

Phenomenologically, more information on the nature of the series can be gained 
if the t dependence is also involved. The well-known (diffractive) dip-bump struc- 
ture of the differential cross section can be roughly imitated by the Glauber (or 
eikonal) series, although more reffned studies within the dipole Pomeron model 
(DP) (linear behavior in ln(s)) @ show that the relevant series is not just the 
Glauber (eikonal) one. A generalization of the DP model including higher terms 
in ln(s) was considered in We mention these attempts only for the sake of 



completeness, although we stick to the simplest case of t = 0, where there are hopes 
to have some connection with the QCD calculations. 

In the present paper we consider a new parametrization for total cross sections 
based on the contribution of a ffnite series of QCD diagrams with relative weights 
(coefficients) and rapidity gaps to be determined from the data. Each set of the di- 
agrams is "active" in "its zone", i.e. the parameters should be fitted in each energy 
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interval separately and the relevant solutions should match. The matching proce- 
dure will be similar to that known for the wave functions in quantum mechanics, 
i.e. we require continuity of the total cross section and of its first derivative. 



2 Description of the model and an example with 
two gluon rungs 

The Pomeron contribution to the total cross section is represented in the form 

N 

^''(^) = E/^(^)^(^-4)^K+'-^)' (1) 

where 

f,is) = j2a,,L\ (2) 

j=0 

So is the prong threshold, 6{x) is the step function and L = ln(s). Here and in the 
following, for s and Sq it is understood GeV^) and So/(l GeV^), respectively. 
The main assumption in Eq. (1) is that the widths of the rapidity gaps In(so) are 
the same along the ladder and are energy independent. Their magnitude is not 



known a priori, but can be related to correlations between jets or multiclusters [p 
formed by single gluons or determined empirically. The functions fi{s) are finite 
polynomials in L, corresponding to finite gluon ladder diagrams in QCD, where each 
power of the logarithm collects all the relevant diagrams. Each time the rapidity 
gap ~ ln(s) exceeds the threshold value In(so), a new prong opens adding a new 
power in L. 

In Eq. (1) the sum over is a finite one, since N is proportional to ln(s), 
where s is the present squared cm. energy. Hence this model is quite different 
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from the usual approach where, in the Umit s oo, the infinite sum of the leading 
logarithmic contributions gives rise to an integral equation for the amplitude. 

To make the idea more clear, we first describe the mechanism in the case of 
three gaps (two rungs) with 

/o(s) = aoo , (3) 

= oio + aii-L , (4) 

f2{s) = a2o + a2iL + 022!/^ ■ (5) 

By imposing the requirement of continuity (of the cross section and of its first 
derivative) one constrains the parameters. E.g., from the equality /i(so) = /o(so) 
the relation 

oio = ooo - On In(so) (6) 
follows. Furthermore, from f2{sl) — fi{sl) one gets 

020 = -2a2i In(so) - 4a22 In^(so) + Oio + 2aii In(so) (7) 

and from the continuity of the relevant derivatives 

021 = dii - 4a22 In(so) (8) 

follows. 

To remedy the effect of the opening prongs and get a smooth behavior at low 
energies, we have included also a Pomeron daughter, behaving like ~ to the 
Eqs. (3) and (4) with parameters 60 and bi respectively (otherwise the continuity 
condition could not be applied to the first, constant term, whose derivative van- 
ishes) . 
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In fitting the model to the data, we rely mainly on pp data that extend to the 
highest (accelerator) energies, to which the Pomeron is particularly sensitive. To 
increase the confidence level, pp data were included in the fit as well. To keep 
the number of free parameters as small as possible and following the successful 



phenomenological approach of Donnachie and Landshoff [|T2|, a single "effective" 
Reggeon trajectory with intercept a (0) will account for non-leading contributions, 
thus leading to the following form for the total cross section: 

at{s)=a''{s) + R{s) , (9) 

where (J^{s) is given by Eq. ([^) and R{s) = as"'^^^^^ (the parameter a is different 
for pp and pp and is considered as an additional free parameter). 

Ideally, one would let free the width of the gap Sq and consequently the number 
of gluon rungs (highest power of L). Although possible, technically this is very 
difficult. Therefore we proceed by trial and error, i.e. make fits for fixed (two and 
three) number of rungs (power of the logarithms). Even within this approximation 
there is some room, as we shall see, to study the sq dependence of the results. 

Notice that the values of the parameters depend on the energy range of the fitting 
procedure. For example, the values of the parameters in /q if fitted in "its" range, 
i.e. for s < Sq, will get modified with the higher energy data and correspondingly 
higher order diagrams included. 

Fits to the pp and pp data were performed from ^/s = 4 up to the highest energy 
Tevatron data (for pp), including all the results from there |T^. To cover the energy 



range with equal rapidity gaps uniformly, sq was chosen to be equal to 144. 

The resulting fit is shown in Fig. 2. The values of the fitted parameters are 



quoted in Table 1. 



3 Three gluon rungs and fits to the pp and pp data 

We cover the energy span available in the accelerator region by four gaps resulting 
in three gluon rungs and consequently as the maximal power. The individual 
Pomeron terms and relevant gaps now are 

4 > < : /o = aoo + &o/s , (10) 

Vso>Vs<{VsQf : fi^aio + bi/s + auL , (11) 

iV^of >V^< iV^of ■ f2 = 020 + a2lL + 022^' , (12) 

(v^o)^ > < (v^o)'' : f3 ^ aso + asiL + a32L^ + a33L\ (13) 

As in the case of three gaps, the individual Pomeron terms and their derivatives 
were matched at the prong values. E.g. oio was determined from the condition 
fo{so) — fi{so), while from the equality of the relevant derivatives, bi can be ex- 
pressed in terms of the other parameters. Ultimately, we are left with nine free 
parameters: ooo, bo, an, 0,22, 032 and 033, each determined in its range, while a and 
a (0) are fitted in the whole range of the data. The parameter sq in principle is also 
free, but as discussed above, we determine it by trial and error, starting with the 
value So — 64. The gap width was chosen such as to cover the whole rapidity span 
by at most four gaps (to have L^). The final value for sq turned out to be 42.5, 
resulting in a sequence of energy intervals ending at y/s — 1800. 

Fig. 3 shows our fit to the pp and pp total cross section data. The values of the 
fitted parameters are quoted in Table 2. The value of the x^/d.o.f. is ~ 1.3, much 
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better than in the case of two gluon rungs. 

The value of the effective Reggeon intercept remains rather low, close to 0.5 
(compare for example with Ref. [0), however interestingly enough its value is 
correlated with the gap width: the smaller the gap width, the higher the Reggeon 
intercept. 

4 Conclusions 

Although high quality fits were not the primary goal of the present study, we may 
conclude that our results are comparable with those of similar analyses |T2|. There is 
still room for some technical improvements in this direction. Our main goal instead 
was to seek for a correct form of the "perturbative" series of total cross sections 
and for regularities in the behavior of the parameters. In fact we find that the 
coefficients in front of leading logarithms in the Pomeron contribution are related 
roughly by a factor 1/10. Notice the alternating signs in front of the logarithms. 
They may reflect the fact discussed in the introduction, namely that each power of 
the logarithms collects various contributions of the same order but from different 
diagrams (see Fig. 1). 

"Footprints" of the prongs at low energies are slightly visible in Fig. 2 (especially 
in the case of pp scattering where the contribution from secondary Reggeons is 
smaller than in pp). A more detailed study of this phenomenon could answer 
the question whether this is an artifact or a manifestation of the Pomeron's basic 
properties. 

The present model and its experimental verification may shed light also on the 
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energy range of the applicability of various approximations to the Pomeron. The 
simplest, Low-Nussinov model with constant cross sections is a crude approxima- 
tion to reality. The inclusion of one gluon rung may be associated with the dipole 
Pomeron. This model has many attractive features, such as selfconsistency with re- 
spect to s-channel unitarity. Note that ln(s) is the strongest rise within the Regge 
pole model. The next order, In^(s), conflicts with the unitarity bound, requir- 
ing that the rise of the cross section does not exceed that of the slope parameter 
(shrinkage of the cone), that in the Regge approach is at most ln(s) (unless special 
assumptions are involved) . This regime seems typical of the Tevatron energy region. 
The role and weight of higher order terms is interesting, but needs more care for 
two reasons: first, too many free parameters make their determination difficult and 
secondly, they violate the Proissart bound, therefore eikonal corrections - otherwise 
present everywhere - here become crucial. A generalization of the above procedure 
within the eikonal model is possible, although the calculations (matching, fitting) 
become more complicated. 

Extrapolations to still higher energies are of great interest. By fitting the model 
to the cosmic ray data and/or future (RHIC, LHC) accelerator data, one could 
explore the role of the new thresholds. On the other hand, from the present fits 
and, hopefully, from the QCD calculations one may try to deduce recursion relations 
and try to extrapolate the value of the coefficients in front of the logarithms. In 
any case, the higher the energy, the more important the unitarity corrections will 
become. Future fits of the model to various data may settle some details left open 
by this paper. 
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TABLES 



So 


bo 


bi 


«(0) 






144 


-74.4 


421 


0.421 


127 


45.9 




«10 


an 


0,20 


021 




35.7 


15.2 


3.44 


57.3 


-5.03 


0.427 



Table 1: Fitted parameters in the case of two rungs. The parameters hi and a... are 
given in units of 1 mb. 





ho 




n(0) 




"PP 


42.5 


-89.3 


43.0 


0.550 


105 


56.3 






On 


O20 


021 


O22 


30.6 


15.5 


3.22 


27.3 


0.0629 


0.211 


030 


031 


032 


033 






106. 


-18.7 


1.68 


-0.0376 







Table 2: Fitted parameters in the case of three rungs. The parameters hi and a... 
are given in units of 1 mb. The value of the x^/d.o.f. is ~ 1.3. 
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Figure 1: Schematic representation of the total cross section in the leading ln(s) 
approximation (first row). Double lines represent protons or anti-protons, vertical 
zig-zag lines are Reggeized gluons, horizontal wavy lines are gluons. The effective 
vertex for two Reggeized gluons and one gluon is defined in the second row. Here 
external lines can represent quarks or gluons. 
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Figure 2: Total cross section calculated up to two gluon rungs and fitted to the 
and pp data. 
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Figure 3: Total cross section calculated up to three gluon rungs and fitted to the 
pp and pp data. 
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